Laser ion acceleration using a solid target coupled with a low density layer 
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We investigate by particle-in-cell simulations in two and three dimensions the laser-plasma inter- 
action and the proton acceleration in multilayer targets where a low density ("near-critical") layer 
of a few micron thickness is added on the illuminated side of a thin, high density layer. This target 
design can be obtained by depositing a "foam" layer on a thin metallic foil. The presence of the 
near-critical plasma strongly increases both the conversion efficiency and the energy of electrons 
and leads to enhanced acceleration of proton from a rear side layer via the Target Normal Sheath 
Acceleration mechanism. The electrons of the foam are strongly accelerated in the forward direction 
and propagate on the rear side of the target building up a high electric field with a relatively flat 
longitudinal profile. In these conditions the maximum proton energy is up to three times higher 
than in the case of the bare solid target. 

PACS numbers: 52.38.Kd, 52.65.-y, 52.65. Rr 

Keywords: PIC simulation; Laser-plasma interaction; Laser-driven ion acceleration 



I. INTRODUCTION 

The availability of laser systems with short and high 
intensity pulses allowed the development of different tech- 
niques for the acceleration of both electrons [l|, Q and 
ions [1, Q exploiting the laser-plasma interaction and the 
high electric field generated in the plasma. The plasma is 
usually obtained from the ionisation of the chosen target 
by the laser pulse itself and, depending on the material, 
different values of the plasma density can be achieved 
leading to different interaction regimes. In the interac- 
tion of an electromagnetic (EM) wave of frequency uj (i.e. 
a laser pulse) with a plasma, the values of the electron 
density rie and of the parameter ric = mea;^/47re^ (so- 
called critical or cut-off density) discriminate between 
two different regimes. For underdense plasmas where 
Ue < ric, according to linear theory the laser pulse prop- 
agates into the plasma and thus the interaction occurs 
through the entire plasma volume. CoUisionless absorp- 
tion of the laser energy may occur, e.g., via the exci- 
tation of plasma waves due to wakcficld generation or 
parametric processes (Raman scattering) 01 • In turn, the 
plasma waves may break and lead to the generation of 
fast electrons, whose energy and number increase with 
the plasma density [5|. For overdense plasmas where 
ne > ric! the laser only penetrates in the "skin" layer 
of thickness ~ c/ujp = (A/27r)-yAi^r/n^ and thus a sur- 
face interaction, rather than a volume interaction occurs. 
A considerable amount of the laser energy is usually re- 
flected, but a sizable or even major part may be absorbed 
via resonant or non-resonant excitation of surface plasma 
oscillations which, in a steep profile, may also break and 
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generate fast electrons [l[ . In this regime usually the ab- 
sorption fraction increases with decreasing density as the 
laser pulse penetrates more deeply into the plasma. 

In the context of experiments using high-intensity 
lasers with optical or near-infrared wavelengths (A = 
0.8 — IfJ-To), most of the times the underdense and over- 
dense regimes correspond in practice to the use of gas 
targets with densities <C ric and of solid targets with 
rig ^ ric, respectively. A regime which is at the boundary 
between underdense and overdense plasma may lead to 
efhcient absorption and fast electron generation, accord- 
ing to the observed scaling with density. Despite the in- 
terest in such "intermediate" conditions with — ric [E~ 
[1| , these conditions have not been accurately investigated 
with experiments because such targets are not straight- 
forward to produce. Actually, in real experiments laser 
prepulses often lead to early plasma formation from solid 
targets so that the short-pulse, high intensity interaction 
actually occurs with an inhomogeneous plasma with both 
underdense and overdense regions. However, such condi- 
tions are usually out of a complete experimental control 
and both interpreting experimental data and optimising 
energy absorption are not straightforward. This issue 
is of particular relevance for proton acceleration via the 
Target Normal Sheath Acceleration (TNSA) mechanism, 
where protons on the rear side of the target are acceler- 
ated by the space-charge electric field generated by fast 
electrons having crossed the target and escaping in vac- 
uum. Thus, increasing the conversion efficiency and the 
energy of electrons generated at the front side yields a 
more efficient TNSA of more energetic protons as well. 

Some recent experiments aimed at improving ion ac- 
celeration by laser pulses have considered solid targets 
of thickness short enough that the long laser prepulse 
produces a near-critical plasma or even ultrashort 
thickness such that the expansion during the interaction 
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leads to self-induced transparency, i.e. to an overdense to 
underdense transition Rather thick (> 10^ /xm) 

foam targets with low density have been used in experi- 
ments either at relatively moderate (~ 10^^ W cm~^) [l2| 
or at very high intensity lO^^W cm~^) [H, and 
relatively "long" (> 5 x lO^fs) pulse durations. Thinner 
foams and ultrashort pulse durations (tens of fs) have 
been studied with particlc-in-ccU (PIC) simulations both 
for intensities in the moderate 10^^ W cm~^) [l^ and 
and uhra-high (~ lO^^W cm"^) ^ regimes. Most of the 
above mentioned studied envision ion acceleration mech- 
anism which may be considerably different from TNSA. 
However, TNSA is the mostly investigated mechanism 
so far, and its peculiar features such as broad spectrum 
and ultra-low emittance proton emission make it most 
suitable for specific important applications such as time- 
resolved proton radiography [l|, |lg . It is thus of interest 
to couple the potential of low-density targets with the 
TNSA scheme. 

The typical proton acceleration scheme in the TNSA 
regime consists of a thin solid foil (Z ~ 1 — lOyum) irradi- 
ated by a high intensity laser (/ = 10^^ — lO^^W/cm ). 
The maximum proton energy is mainly limited by the 
laser energy absorbed by the electrons. Material science 
offers opportunities for manufacturing targets with low 
and controlled density, such as foam targets. The pos- 
sibility to couple a thin solid foil with a near critical 
density layer attached on the irradiated side can enhance 
the laser energy absorbed by the target with a consequent 
increase of the electron energy and the accelerating field 
arising on the rear side of the target. An "advanced" 
TNSA regime may be obtained with a higher proton en- 
ergy considering the same laser characteristics. 

Here wc report a simulation study of laser absorption 
and ion acceleration for a three-layer target configura- 
tion, where a low density, few-microns thick ("foam") 
layer and an ultrathin ( "impurity" ) proton layer are de- 
posited on the front and rear sides, respectively, of a 
("solid") foil target. A similar three- layer target config- 
uration has been studied by Nakamura et al. [Tsj using 
2D simulations at normal incidence. An enhancement of 
proton acceleration due to the presence of the foam has 
been shown together with the effects of the field ionisa- 
tion. In the present work we extend the investigation in 
different directions. We present fully three dimensional 
(3D) simulations and 2D simulations with oblique inci- 
dence. We use the PIC code ALaDyn [13, 111, in both 
2D and 3D, to investigate the dynamics of the laser inter- 
action with the slightly overcritical density plasma and 
the role of the foam electrons in the rise of the longitu- 
dinal electric field which accelerates the protons. Both 
3D and 2D simulations show the proton energy "gain" 
obtained considering a target with a foam layer instead 
of a "bare" target with no foam. We show the impor- 
tance of a 3D study to correctly determine the maximum 
value of the proton energy and the dynamics of the accel- 
eration. A wide paramctrical study in 2D as a function 
of foam thickness and density and of the incidence angle 



of the laser is also reported. We quantitatively analyze 
the absorption of the laser energy by the foam-solid sys- 
tem, considering different parameters of the target and 
relating it to the energy of the accelerated protons. We 
address in detail an important aspect of the proton accel- 
eration mechanism, namely the formation of a strongly 
non-equilibrium electric field at the rear side of the tar- 
get, due to the prompt escape in vacuum of bunches of 
highly relativistic electrons created in the laser foam in- 
teraction. 



II. SIMULATION SET-UP 

In the PIC simulations we consider the case of a p- 
polarised laser pulse incident (at normal incidence) on 
an already ionised, ideal plasma. The target is composed 
by three layers: a thin metallic foil (~ 0.5/im) (main 
layer), a thicker low density foam (1 — 12^m) attached 
on the side directly irradiated by the laser beam, and 
a third layer representing the thin (^ 50nm) contami- 
nants layer on the rear surface. The ionization state of 
the target layers is fixed and the charge over mass ratio 
arc Z/A = 1/2, Z/A = 1/3 and Z/A = 1 for the foam 
layer (e.g. C^+), metal foil (e.g. Al®^) and contami- 
nants (H^) respectively. The laser propagates along the 
X direction and in every case considered independently 
on the foam thickness, it reaches the maximum focusing 
on the front surface of the metal foil. In all the simu- 
lations reported in the paper the laser pulse wavelength 
was fixed to A = 0.8/im and the pulse shape was kept as 
sin^ and Gaussian in the longitudinal and the transverse 
direction respectively with t = 25fs and wq = Sfim being 
the time duration (FWHM) and pulse waist. The peak 
value of the dimensionless amplitude of the laser pulse 
ao — 8.5-10^^°A[yum]/^/-^[W/cm^] was varied in the range 
ao = 3 — 20; which corresponds to a peak power range 
P = 2.8 - 128TW and total energy U = 0.075 - 3.2J, 
well within the capabilities of many present-day facil- 
ities. The spatial resolution is Ax = c/2u!p ~ A/111 
Ay = Az = c/ujp, where Wp is evaluated from the maxi- 
mum density of the plasma (i.e. the solid foil density). 

Two 3D simulations have been performed considering 
the same solid target with and without a foam layer 
Uf = 271c, If = 2yum. In these cases the "contami- 
nant layer" and metallic foil density and thickness are 
Ur = 9nc, Ir — 0.05/im and Um = 40nc, Im — 0.5/im re- 
spectively. The electron population of the main foil and 
of the foam layer have been sampled with respectively 36 
and 8 macro-particles/cell, whereas the ion populations 
with 12 and 8 macro-particle/cell; protons and electrons 
of the contaminants layer are sampled with 36 macro- 
particles/cell. The total number of grid points of the 
3D simulations has been about 3 • 10^ and the number 
of macro-particles is about 6 • 10^; these simulations re- 
quired about 30000 CPU-hours on the SP6 machine of 
CINECA (Bologna, Italy). 

A more extensive investigation spanning on a wider 
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range of foam parameters has been performed with sev- 
eral 2D simulations. The electron density of the foam 
Tif and its thickness If have been varied: = 1 — Stic, 
If = 1 — 12/xm, whereas the main layer density has been 
fixed to rim = SOric. The nmnber of macro-electrons per 
cell was 81/cell in the solid density layer, 9 for the con- 
taminants layer and 16 or 25 for the foam layer, whereas 
the number of macro- ions of the three species was 9/cell. 

The time t = corresponds to the instant when the 
laser starts to interact with the plasma; most of the sim- 
ulations have been stopped after 50/^m/c = 166 fs. A 
comparison of the results obtained in 2D and 3D shows 
how in 3D the acceleration process at i = 166 fs is almost 
finished, whereas in the 2D simulations the acceleration 
process is still significant even after nearly 400fs. In a 
2D simulations the "point" charges effectively behave as 
infinite wires with finite linear charge density. On longer 
time-scales the problem is nearly pure electrostatic and 
the logarithmic behaviour of the potential ~ — ln(r)) 
leads to a much slower "saturation" of the acceleration 
process and higher proton energies. 



III. RESULTS AND DISCUSSION 

A. 3D simulations 

The density of the solid foil chosen for the 3D simu- 
lation is Tim ~ 40ric; though less than the value for a 
real solid metal (several hundreds ric), it is considerably 
higher than the critical density and, for a peak inten- 
sity ao = 10, it ensures that the plasma is overdense 
even accounting for relativistic transparency effects [l^ . 
A set of 2D simulations strengthened the choice of the 
foil density: we considered cases varying the laser in- 
tensity (flo = 10 — 20) and the plasma density of the 
main foil (n^ = 40 and SOric). For the lower laser in- 
tensity (ao = 10), when no foam layer is present the 
results obtained with the two different main target den- 
sities {rim = 40 or SOric) are similar, (proton maximum 
energy is 11 instead of 9 MeV) but a lower density, as 
expected, leads to better absorption, whereas they are 
barely distinguishable in presence of a foam. For the 
higher value of the pulse intensity (ao = 20), the dif- 
ference between the case of rim = 40nc or rim — SOric is 
more significant with the lower density leading to consid- 
erably higher proton energy and enhanced laser energy 
absorption in both cases, with and without foam. This 
can be explained by the fact that for the high intensity 
case the foam layer considered is too thin to absorbs a 
sizable fraction of the total laser energy, thus the interac- 
tion mostly occurs with the solid foil and is more sensitive 
to the density of the latter. Additional simulations have 
been performed doubling the thickness of the solid foil 
{rim = SOric, Im = l/^rn instead of Im = 0.5/im). The 
maximum proton energy in presence of the foam layer 
{If = 4/zm rif = Iric) is reduced by a few percent (~ 7%, 
24 instead of 26 MeV), whereas for a bare target, as also 
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FIG. 1. a): proton energy spectrum at t = 166 fs considering 
a laser pulse with ao = 10 and a target with a foam layer 
71/ = 2nc, If — 2/i. b): proton maximum energy evolution 
with respect to time: comparison 3D (solid) and 2D (dashed) 
cases without foam (red) and with (red) n/ = 2nc, If = 2/i, 
n,-n = 40nc, ao = 10. 



reported in previous work [20j , it decreases more substan- 
tially (^ 25%, 6.5 instead of 9 MeV). It can be concluded 
that if the foam thickness is sufficient to absorb a major 
part of the laser energy, the characteristics of the solid 
foil are less crucial. 

The simulations led to proton energy spectra with an 
exponential profile and a cut-off typical of the TNSA 
regime (Figure [T]). The energy spectra arc obtained con- 
sidering all the protons of the contaminants layer achiev- 
ing a spatial integration in the transverse direction. The 
presence of a collimator after the target would then se- 
lect the high energy tail of the spectrum similarly to 
the TNSA case. At a laser intensity corresponding to 
ao = 10, in presence of a foam layer, although rather 
thin (2/im), the maximum proton energy Emax is much 
higher, Emax.f — 14MeV than without foam Emax,b — 
6McV (see fig[T]). The same simulations have been per- 
formed also in 2D: the cut-off energy of the protons is 
overestimated (by a factor of about 2 at time 166fs) but 

the ratio E^Zxjl E^rZx,, ^ E^Zxj I E^ax,b ^ 2.3 is pre- 
served. It is therefore evident that the 3D analysis is 
essential to evaluate the maximum ion energy quantita- 
tively. Nevertheless, the observed dynamics of the laser- 
foam interaction is qualitatively similar in 2D simula- 
tions, thus we use the latter for a parametric study aimed 
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at both showing the existence of an optimal foam thick- 
ness as a function of the laser amplitude and angle of in- 
cidence, and at evidencing features of the electron accel- 
eration and sheath field formation processes (Sec lIII B|) . 

In Figure [2] the volume rendering of the electron den- 
sity is presented together with the distribution of the 
longitudinal electric field. The electrons from the foam 
arc accelerated in the forward direction and reach the 
rear side of the target. The electron cloud expands with 
a roughly spherical symmetry for several microns and 
displays regular structures in the longitudinal direction 
which are more evident in the 2D plots and will be fur- 
ther discussed. The resulting charge separation leads to 
a strong longitudinal electric field which exceeds 5 TV/m 
and extends in the longitudinal direction according to the 
electron distribution for about 10 microns. The slice dis- 
played in the Figured] highlights the longitudinal electric 
field of the laser pulse and shows how the laser reaches 
the solid foil and is then reficcted by the high density 
plasma. Thus, even if rie > ric the foam plasma is effec- 
tively underdense for the laser pulse, due to relativistic 
penetration and ponderomotive channeling effects. Since 
an extended analysis of the laser plasma interaction and 
ion acceleration for several values of the parameters can 
be hardly performed in 3D we continued the investigation 
with 2D simulations. 



B. 2D simulations 

The detailed analysis of the results for several target 
parameters has been carried on the basis of 2D simula- 
tions and the discussion will now focus on 2D cases only. 
In the light of the 3D results and the comparison with 
the analogous 2D cases, we will focus the attention on the 
main features of the laser-plasma interaction as well of 
the foam target parameter rather than on a quantitative 
estimation of the proton energies. 

When a low density foam layer is present, the en- 
ergy absorption mechanism is different from the case of a 
highly overdcnsc plasma: in the range of parameters here 
considered the laser propagates through the foam and is 
not effectively reflected until it reaches the solid layer. A 
minor part of laser energy is then absorbed by the high 
density plasma similarly to the case without foam. Dur- 
ing the interaction with the foam, the laser pulse accel- 
erates the electrons to relativistic velocities and its en- 
ergy is considerably depleted. We analysed the energy 
balance of the EM energy and kinetic energy of the par- 
ticles (Figure [3]). The parameter scan in 2D, Figure [3]-a, 
shows that with the increase of the foam thickness (If), 
for a constant foam density, the ratio of the energy ab- 
sorbed by the electrons over the initial laser pulse energy 
increases. For a constant foam density nf = ric, the es- 
timated reflectivity is reduced from 92% for the case of 
"bare" solid target, to 18% using a 8^m foam. The right 
panel of Figure [3] shows the time evolution of the particles 
kinetic energy for a case with a foam layer with thickness 




FIG. 2. 3D simulations (ao — 10, wq = 3/im, If = 2/im, 
rif — 2nc): electron density in logarithmic colour scale (a) 
and longitudinal electric field (b) at t = 66 fs in presence 
of the foam. Half of the simulation space has been removed 
(2 > 0) for easier viewing. 



If = 8/im and density n/ = Uc- When a sufficiently thick 
foam is present {If > 4/zm), the electrons gain a notably 
large fraction of the initial laser energy (up to > 50% of 
the total energy whereas without foam the correspond- 
ing value is about 5%). The ions slowly gain energy at 
the expense of both the electron kinetic energy and, to 
a lower extent, the electrostatic energy. It is remarkable 
how the protons from the contaminants, although form- 
ing a thin (50 nm) and low density (n^ = 9nc) layer, get 
about 50% of the total kinetic energy absorbed by all 
the ion population of the simulation, accounting for up 
to 10% of the initial laser energy (to be compared with 
~ 1% of the corresponding case without foam). 

The cut-off value of the proton energy spectra (see Fig. 
2]) is sharp and is related to the maximum value of the 
field acting on the protons at the beginning of the accel- 
eration process, analogously to the ordinary TNSA case 
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FIG. 3. a): maximum value of the electron energy ab- 
sorbed by the target for different cases, b): for the case with 
I f = 8^m: the time evolution of the energy normalised to the 
initial laser energy of: all electrons (black), foam electrons 
only (red), all ions (green), contaminant protons only (blue). 



[21| . A summary of the maximum proton energy ob- 
tained at time t ~ 166 fs in the different cases analysed 
by the 2D simulations is reported in Figure [S]-a as func- 
tion of the foam areal density (thickness in /.im times 
Uf /hc). This comparison shows how at a given laser in- 
tensity, there is an "optimal" foam thickness for each 
value of the density (for oq = 10, 8 — 12/im for ti/ = n-c, 
2 — 4^m for n/ = 2nc and 1 — 2^m for n/ = 4nc) and, at 
least in the range of parameters considered, this corre- 
sponds to the same value of the areal density. In Figure 
^}p) the dependence of the maximum proton energy with 
respect to the foam density is reported considering laser 
pulses of different powers. The plot shows that at a given 
laser intensity, an optimal value of the foam areal density 
can be found and it is approximately proportional to the 
laser intensity. 

The increased proton energy is a direct consequence of 
the stronger electric field arising if a foam layer is present. 
The accelerating field at the rear side of the target is 
generated by the charge displacement of the electrons 
which are accelerated by the laser and propagate in the 
forward direction. After the laser pulse has been reflected 
by the solid layer {t > 70fs), in presence of a foam, the 
longitudinal electric field exhibits a different shape and 
a maximum value about 3 times higher if compared to 



FIG. 4. a): proton energy spectra ed t = 166 fs/ considering 
a laser with ao — 10 and foam layer with density n/ — ric- b): 
maximum proton energy evolution with respect with time for 
different foam thicknesses (ao = 10, n/ = ric). 



the case of a bare solid foil. Figure [6] shows the line- 
out of the longitudinal electric field on axis for the cases 
with foam (If = Sfim Uf = ric) and without for ag = 
10. An exponential decrease is accompanied by a nearly 
uniform field for several microns from the rear surface. 
Whereas the exponential decrease can be attributed to 
a "hot" electron population which expands around the 
target similarly to the ordinary TNSA ease. In this case, 
the electric field is generated due to the highly relativistie 
electrons promptly escaping far away from the target, 
and it is strongly different from the expression that is 
obtained in the assumption of a Boltzmann equilibrium. 
For foam thicknesses lower than the case of Figure [6] this 
"step-like" structure of the electric field is still present 
but not as clearly distinguishable and more similar to 
the case of "ordinary" TNSA fields. 

The electron spectra shown in FigiT] for ao = 10 and 
foam densities Uf = 2nc, Uf = ric are characterized 
by cut-off energies of ~ 30, ~ 50 MeV and rough es- 
timates of the electrons "temperature" of ~ 9 MeV 
and ~ lO.SMeV, respectively. These values are much 
higher than the widely used "pondero motive" sc aling of 
the electron temperature Tp = meC^( \/l + a^/2 — 1) ~ 
6.1meC^ =3.1 MeV for interactions with highly over- 
dense plasmas, where the fast electrons have energies of 
the order of the "quiver" energy in vacuum, and the "ac- 
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celeration length" is of the order of a wavelength in such 
case. The comparison shows that for the foam targets 
the production of fast electrons occurs along the volume 
of the foam and suggests that peculiar mechanisms of 
electron acceleration are at play. In Figures [8] the struc- 
tures of the expanding electrons show an antisymmetry 
with respect to the propagation axis. The antisymmetry 
is related to the polarisation of the pulse, as it appears 
for a p-polarisation only (i.e. for E in the simulation 
plane) An analogous simulation has been run consider- 
ing the same conditions but with an s— polarised laser 
pulse instead. The "bunching" of the electrons is still 
present but less evident and the electron density is per- 
fectly symmetric with respect to the propagation axis. 
Moreover, the electrons reach much lower energies with 
and energy spectrum which extends up to about 16 MeV 
instead of 32 MeV with a p— polarised laser pulse. As a 
consequence also the maximum proton energy is lower: 
17 MeV instead of 30 MeV. The comparison gives evi- 
dence of a strong contribution of the transverse electric 
field to electron acceleration. Again, eSicient penetra- 
tion into the overdense plasma is necessary in order for 
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FIG. 5. a): Proton maximum energy for ao — 10, different 
foam tfiicknesses and foam densities (black squares: Uf = Uc, 
red circles: n/ = 2nc, blue triangles: n/ = 4nc) as function of 
the areal density {IfUf /ric). b): Proton maximum energy ver- 
sus foam density for different thickness (2 ^m solid line, 4 fim 
dashed line) and laser intensity (ao = 3 magenta pentagons, 
ao = 7 green diamonds, ao = 10 black squares, ao = 14 red 
circles, oo = 20 blue triangles). 



field (t=100fs) 



5 





1 lf=8, n,=nj! 




1 no foam 










w 





-5 5 10 15 

X [urn] 



FIG. 6. Cut of the longitudinal electric field on the focal axis 
z = with a foam layer(?i/ = In^ If = 8fim) and without, for 
ao = 10 at f = 100 fs after the interaction's onset. The peak 
values are Ex,m,ax = 5.4TV/m and Ex.max = 1.8TV/m. 

the coupling with the transverse field to be effective at 
normal pulse incidence. In the case of a highly overdense 
plasma, such that the laser penetrates only in the skin 
layer, at normal incidence the coupling with the trans- 
verse electric field is inefficient and absorption is domi- 
nated by "J X B heating" 

The evidence of electron acceleration over the volume 
of the foam, bunching of electrons and more efficient 
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FIG. 7. Electron spectra. Electrons from the foam layer 
(uf = ric a), 71/ = 271c b) and electrons coming from the metal 
foil (c) for different simulations using a foam layer 71/ = 2nc- 
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FIG. 8. P— polarised laser pulse (ao — 10) incident on a target 
with foam layer nf = luc, If = S/im at normal incidence: 
electron density (a) and longitudinal electric field (b). 



to shorter pulses for reasons of numerical feasibility and 
relevance to the parameters of already proposed experi- 
ments, however it appears of interest for further study to 
test the scaling of the foam-enhanced acceleration mech- 
anism at higher energies of the laser pulse. In this per- 
spective, from the technical point of target manufactur- 
ing the deposition of a foam layer may be easier than 
engineering a hollow microcone and also relax alignment 
and pointing requirements. 



C. Oblique incidence 

In most of the TNSA experiments, the laser pulse is fo- 
cused on target at oblique incidence with p— polarisation. 
This increases the absorption and usually leads to higher 
proton energies. When a foam layer is present on the 
irradiated side of the target, the proton energy is less 
sensitive to the incidence angle. We performed some 
2D simulations considering targets with a rather thin 
foam layer (If ~ 2/im) of different densities and a 
polarized laser pulse incident on target at different an- 
gles (0°, 30°, 45°, 60°, Wo = S^im, ao = 10, r = 25fs). 
In the case of bare target the results show how, despite 
the effective intensity "on target" is reduced by a fac- 
tor 1/V2 (for 45°), the elect ron heating mechanisms at 
oblique incidence is more efficient [l|, [2^ leading to higher 



coupling for P-polarization, combined with the observed 
penetration and channeling of the laser pulse inside the 
foam, suggests that the dominant mechanism of fast elec- 
tron generation may be similar to that observed in hol- 
low microcone (funnel-like) targets, which allowed to ob- 
tain the highest proton cut-off energy experimentally ob- 
served to date |22l | . The mechanism, named "direct laser- 
light-pressure acceleration", relies on the effective local 
grazing incidence of the laser pulse on the microcone 
walls, where the P-component of the electric field ex- 
tracts electrons which are then accelerated by the com- 
bined action of ponderomotive force and self-generated 
fields, resulting in electron temperature much higher than 
the ponderomotive scaling, the increase observed in 2D 
simulations of Ref. [2^ at ao = 12 from Tp = 3.9 meV to 
9.2 meV being similar to what we observe for the foam 
target. In this latter case, the penetration of the laser 
pulse inside the foam effectively yields a self-generated 
microcone or funnel, providing a similar coupling at graz- 
ing incidence with the channel walls. The comparison 
with experiments and simulations in Ref. (22j may be only 
qualitative because of both the prepared microcone struc- 
ture and the considerably longer and wider laser pulse. 
This implies that the total pulse energy is higher by some 
two orders of magnitude than our case, which may ac- 
count for the higher proton energy at similar values of 
the pulse amplitude. In the present paper we restrict 




FIG. 9. Longitudinal phase space of the electrons of the foam 
(n/ = Iric, If = 8/im, ao = 10) at two different times; px in 
normalised rricC units. 
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FIG. 10. P— polarised laser pulse (ao = 10) incident on a 
target with foam layer n/ = 2nc = 2fim with a 45 ° an- 
gle of incidence at time t — 83fs: electron density (a) and 
longitudinal electric field (b). 



proton energy (12 vs. 9 MeV). In presenee of a foam the 
proton buneh obtained is not anymore symmetric with 
respect to the target normal and its propagation direc- 
tion is slightly tilted (~ 2°) off axis. One simulation 
has been run considering a 45 ° incidence and a thick 
foam {If = 8fj.m and rif = ric), which was an optimal 
case for normal incidence, and the maximum proton en- 
ergy decreased to less than 20 MeV compared to the 30 
MeV of the case at degrees. On the other hand when 
a thinner foam is considered If ~ 2fim and rif = 2nc 
(Figure llip the laser-plasma coupling is more effective 
and the resulting maximum proton energy is comparable 
to the ease of normal incidence. Whereas for an angle 
of incidence of 30°, the results are quite similar to the 
ease of normal incidence, although the trend of the max- 
imum proton energy differs for higher values of the foam 
density, increasing the angle of incidence, the maximum 
proton energy become very sensitive to variations of the 
foam density and the optimal value is found in a narrow 
range of values see Figure [TT] 



IV. CONCLUSIONS 

We presented an analysis of a laser driven ion accelera- 
tion mechanism considering a target configuration where 
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FIG. 11. Survey of the maximum proton energy at t — 166fs 
considering different target densities and angles of incidence. 



a solid density thin plasma is coupled with a low density 
layer attached on the irradiated side. We extended the 
study of Nakamura et al. in several directions: per- 
forming 3D simulations, 2D simulations on a wide range 
of parameters and studying the effects of oblique inci- 
dence; moreover we discussed the mechanism of electron 
acceleration and found analogies with the work of Gail- 
lard et al. [22|. The presence of the low density layer 
strongly increases the energy absorption of the target, 
compared to the case of a bare solid foil. For a laser with 
a peak power of 32TW, the energy absorbed by the elec- 
trons goes from about 5% of the total laser energy for 
the case without foam to over 60% for the cases where a 
foam layer with a thickness of about Afim is considered. 
This translates in a much higher number of energetic elec- 
trons with a considerably higher energy; if a foam layer is 
present the electron maximum energy is more than three 
times higher and the number of fast electrons is up to 
one order of magnitude higher if compared to the case of 
a bare solid target. These fast electrons, escaping from 
the rear side of the target, build up a strong electro- 
static field which easily exceeds lOTV/m and display a 
different. During the interaction with the foam, the laser 
propagates through the low density plasma and acceler- 
ates to relativistic velocities the electrons of the foam in 
its direction of propagation, differently than what hap- 
pens in the surface interaction with a highly overdense 
plasma, to relativistic velocities. The comparison of 2D 
simulations with either S- and P-polarization at normal 
incidence shows that electron acceleration is enhanced 
for P-polarization, which allows a coupling of the laser 
electric field with the walls of the laser-drilled channel. 
These features and the relative enhancement of the elec- 
tron temperature with respect to the typical "pondero- 
motive" scaling of laser-solid interactions also suggest the 
acceleration mechanism to be similar to the one occurring 
in hollow microcone targets. 

The fast electrons propagate through the target and 
reach the rear side creating a strong, quasi-uniform elec- 
tric field. This contribution adds to the electrostatic field 
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arising from the "thermal" expansion of the hot electron 
population from the solid foil, which is also irradiated by 
the laser similarly to the "pure" TNSA regime. This 
configuration has been tested by means of PIC simu- 
lations in 3D and 2D. The 3D simulations allowed for 
a more quantitative estimation of the proton accelera- 
tion and showed how even if a relatively thin (If = 2/ini 
Tif = 2nc) foam layer is added on the front surface of a 
sohd target {If =0, 5/xm rif = AQric) the energy absorp- 
tion is strongly increased and the maximum proton en- 
ergy obtained considering a 32TW laser pulse (r = 25fs, 
wq = 3/xm, with a peak intensity of 2.1 10^°W/cm^), 
reaches 15 MeV which is more than two times higher 
than the case without foam layer (6.5 MeV). The 2D sim- 
ulations, on the other hand, allowed for more extensive 
investigation of the regime changing the target parame- 
ters (for different thicknesses // = 1 — 12/im and densities 
Uf = 1 — Sric of the foam layer) and the laser peak in- 
tensity over several values. For each laser power, ranging 
from 3 up to 128 TW, an optimal foam thickness can be 
found at a given density and the maximum proton energy 
can be several times higher than case of bare solid target. 
The presence of the foam proved to strongly reduce the 
role played by the areal density of the solid target for the 
proton acceleration. To maximize the efBciency of the 
process, it is not necessary to reduce the thickness of the 



solid target to very low values (<C 1/im) hence avoiding 
the need to use targets prone to be damaged by the pre- 
pulse. If the laser irradiates the target at moderate angles 
(30 °) the proton energy is not strongly affected, whereas 
for larger angles the oblique incidence is considerably less 
efficient in producing high energy protons. The picture 
arising from this investigation is richer than the case of 
a pure TNSA regime and an experimental investigation 
using high contrast laser might allow to study the effects 
of near critical plasma with a tighter control if compared 
to "pre-heated" expanding targets. 
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